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ABSTRACT 

Using three-dimensional (3D) simulations of neutrino-powered supernova explosions we show that the hy- 
drodynamical kick scenario proposed by Scheck et al. on the basis of two-dimensional (2D) models can yield 
large neutron star (NS) recoil velocities also in 3D. Although the shock stays relatively spherical, standing 
accretion-shock and convective instabilities lead to a globally asymmetric mass and energy distribution in the 
postshock layer. An anisotropic momentum distribution of the ejecta is built up only after the explosion sets 
in. Total momentum conservation implies the acceleration of the NS on a timescale of 1-3 seconds, mediated 
mainly by long-lasting, asymmetric accretion downdrafts and the anisotropic gravitational pull of large inho- 
mogeneities in the ejecta. In a limited set of 15 M Q models with an explosion energy of about 10 51 erg this 
stochastic mechanism is found to produce kicks from <100kms _1 to >500kms _1 , and >1000 kms -1 seem 
possible. Strong rotational flows around the accreting NS do not develop in our collapsing, non-rotating pro- 
genitors. The NS spins therefore remain low with estimated periods of -500-1000 ms and no alignment with 
the kicks. 

Subject headings: supernovae: general — pulsars: general — stars: neutron — stars: kinematics and dynamics 



1. INTRODUCTION 

Young neutron stars (NSs) possess average space velocities 
around 400kms~ 1 , much larger than those of their progeni- 
tor stars, implying that they are acc elerated during the birth in 
a supernova (SN) explos i on (e.g., i Faucher-Giguere & Kaspil 
120061: iHobbset all 120051: lArzoumanian et all 120021) . More- 
over, alignment of the N S spin and kick was inferred for 
the C rab and Vela pulsars dKaplan et ali 2008: N g & Romanil 
120071) and several other young pulsars (see [Wang et al] |2006l 
and references therein) from comparisons of the direction of 
proper motion with the projected rotation axis as determined 
from the symmetry axis of the pulsar wind nebula on X-ray 
images. The same conclusion was drawn for samples of radio 
pulsars from the linear polarization of t he pulses, whose po- 
sition angle r eflects the spin direction (Joh nston et al.ll2005t 
iRankinl 120071) . However, the Crab and Vela pulsars may not 
be good cases for determining misalignments because they are 
moving at smaller speeds than the average pulsar population 
and therefore the unknown velocity of the progenitor implies 
a bigger uncertainty. On the other hand, also the radio data do 
not seem to make a clear case for a general alignment of spin 
and kic k directions in the ref erence frame of the progenitor's 
motion dJohnston et al .1120071) . 

Analysing the observ ational informati on, also of charac- 
teristics of NS binaries, lLai et al.l (120011) concluded that the 
NSs received their kicks most probably at the time of the 
SN explosion. A large variety of mechanisms for natal 
kicks has been proposed, either by hydrodynamical effects 
linked to large-scale asymmetries of the SN explosion (e.g. , 
Herantlll995[lBurrows & Haveslll996t Uanka & Mulled ll994t 



Thompson 2000; Scheck et al. 2004, 2006) or by anisotropic 
neutrino emission from the nascent NS (e.g., Chugail Il984t 
iBurrows & Wooslevl[T986l ISocrates et al.ll2005l) . However, it 
is very difficult to produce even only a one-percent global 
dipole asymmetry of the neutrino emission, which is needed 
for a kick of 300 km s _I . For this to be possible one has to in- 
voke controversial assumptions like very strong global dipo- 
lar magnetic fields inside the NS (> 10 16 G; e.g.. lArras & Lail 



1999), arguable neutrino properties ( e.g., sterile neutri nos, 
large neutrino magnetic moments; e.g., Fuller et al. 20031). or 
unsettled mechanisms to create str ong emission asymm etries 
in the neutrinospheric region (e.g. . ISocrates et afll2005l). 

On the basis of 2D SN models IScheck et al.l (120041120061) 
argued that the standing accretion shock instability (SASI ; 
iBlondin et al.l2003tlFoglizzo & Taggerl2000HFoghzzol2002h . 
which grows after shock stagnation and causes large global 
non-radial asymmetry of the accretion flow to the NS and of 
the beginning SN explosion, can lead to kicks of typically sev- 
eral hundred kms 1 and even more than 1000 km s" 1 if the 
dipole ( I — 1) component of the asymmetry was sufficiently 
strong. Blondin & Mezzacappj (120071) showed in idealized, 
stationary-accretion setups that SASI spiral modes may also 
have the potential to generate pulsar spin periods consistent 
with observations. In this Letter we present the first simula- 
tions of SN explosions from core bounce to -1.5 s later that 
confirm the potential of SASI-induced asymmetries to pro- 
duce typical NS kicks in the more realistic 3D environment of 
collapsing stellar cores. 

2. NUMERICAL SETUP 

We use the explicit finite-volume, Eulerian, multi- 
fluid hydrodynamic s code Prometheus (iFrvxell et al.l 1 199 it 
Mulle retal] 1 1991 alibi) . The advection of nuclear species 
is treated by the Consistent Mul ti-fluid Advection (CMA) 
scheme of Plewa & Muller (1999), and self-gravity accord- 
ing to Muller & Steinmeta (119951) . General relativistic cor- 
rections of t he m ono pole are taken into account as in 
IScheck et all ({2006) and lArcones et al l (120071) . 

All simulations are computed on an axis -free overlapping 
"Yin- Yang" grid (Kagev ama & Satol 120041) in spherical po- 
lar c oordinates, which was recen tly implemented into our 
code (fWongwathana rat et alJl2010h . We use 400(r) x 47(0) x 
137(0) x 2 grid cells corresponding to an angular resolution 
of 2° and covering the full An solid angle. Better resolution in 
the inner region of our computational domain is achieved by a 
constant radial zone size of about 0.030km up to r a 100km 
The radial grid is logarithmically spaced beyond this radius 
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TABLE 1 

Explosion and NS properties for all models at 1.4 s after bounce. 



Model M ns f cxp 

1 spin 

[M e ] [ms] [B] [km/s] [km/s 2 ] [gcm 2 /s] [°] [ms] 

W15-1 1.37 246 1.12 331 175 1.51 117 652 

W15-2 1.37 248 1.13 405 144 1.56 58 632 

L15-1 1.58 421 1.13 161 66 1.89 148 604 

L15-2 1.51 381 1.74 78 3 1.04 62 1041 
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Fig. 1 . — Time evolution of the SN explosion energy £ cxp (top) and of the 
average shock radius i? s for all computed models. Since models W15-1 and 
W15-2 differ only in the initial random seed perturbations of the radial veloc- 
ity field, the results for these models are essentially identical and the short- 
dashed and solid lines lie on top of each other. Despite the similarity of the 
global parameters, both models develop different explosion asymmetries and 
neutron star kicks (see Table[T]and Fig. [2). 

up to an outer grid boundary of R b = 18000 km, which is 
sufficient to prevent the SN shock from leaving the computa- 
tional domain during the simulated time. At R b a free outflow 
boundary is chosen. 

The high-density inner core of the proto-NS (PNS) is ex- 
cised and replaced by a point mass at the coordinate ori- 
gin. The shrinking of the PNS is mimicked by a retraction 
of the closed inner boundary at Rn, together with the radial 
grid. Hydrostatic equilibrium is assumed at R&. The stel- 
lar fluid is described by the tabulated microph ysical equa- 
tion of state (EoS) of Uanka & Miilleri d 19961) . Neutrino 
transp ort and neutr i no-ma tter interactions are approximated 
as in iScheck et al.1 (120061) by radial integration of the one- 
dimensional (spherical), grey transport equation for all an- 
gular grid directions (8, <p) independently. This "ray-by-ray" 
approach allows for angular variations of the neutrino fluxes. 

3. INVESTIGATED MODELS 

Our models W15 and L15 are based on two non-ro tating 
15 M progenito rs, s!5s7b2 oflWooslev & Weaver! d 19951) and 
a star evolved by Limon gTet al.l ( 120001) . which were followed 
through collapse to 15 ms after bounce with the Prometheus- 
Vertex code in one dimension (R. Buras and A. Marek, pri- 



vate communication). To break spherical symmetry, random 
seed perturbations of 0.1% are imposed on the radial veloc- 
ity (v r ) field. Explosions with chosen energy are initiated by 
neutrino heating depending on suitable values of the neutrino 
luminosities imposed at the lower boundary. 

Our models W15-1 and W15-2 differ only by the initial 
seed perturbations, while L15-2 is set to have a higher ex- 
plosion energy than L15-1. Table Q] lists the corresponding 
explosion energies E exp at ~ 1.4 s p.b. and explosion times f exp 
defined as the instant when £ e xp (i.e., the total energy of all 
zones where the internal plus kinetic plus gravitational energy 
is >0) is 10 48 erg. Figure Q] shows that this moment coincides 
roughly with the time when the average SN shock radius, R s , 
exceeds 500 km. The W15 models explode earlier and their 
NS (baryonic) mass, M ns (defined as enclosed mass at a radius 
R ns where the density is 10 11 gem -3 ), is therefore smaller. 

4. NEUTRON STAR KICKS AND SPINS 

Neutrino heating around the NS triggers delayed explo- 
sions, which are preceded (and supported) by a postbounce 
phase of several 100 ms of violent SASI and convective activ- 
ity. This leads to large-scale asymmetries of the ejecta, po- 
tentially imposing a kick and spin to the NS, for which we 
evaluate our simulations in a post-processing step. The time 
evolution of the NS kick velocity v ns , NS angular momen- 
tum 7 ns , lateral and longitudinal angles of the kick direction 
(#kkk and 0kick, respectively), and the angle t? s k between spin 
and kick directions are plotted in Fig. [2] The final numbers 
are given in Table[T](7 ns values there are normalized by 10 46 ), 
supplemented by the NS acceleration a ns at the end of the runs 
and our estimated NS spin periods r sp i n . 

Since the NS is excised and replaced by a point mass at the 
grid center, it cannot travel. Nevertheless, like a wall reflect- 
ing a bouncing ball, it can absorb momentum (but will not 
move as if it had an infinite inertial mass). Because of total 
linear momentum conservatio n in the rest fram e of the pro- 
genitor, we can (following Schec keit al J 120061) compute v ns 
from the negative of the total momentum P gas of the gas out- 
side of the NS as 

V„s(f) = -Pgas(0/M ns (f), (1) 

where P gas = J R d V pv. Tests confirmed very good linear 
momentum conservation of our code, while angular momen- 
tum is more difficult to conserve, e.g., when a rotating gas 
mass is in rapid motion across large distances on the grid. We 
therefore estimate J ns again as the negative of the angular mo- 
mentum of the exterior gas, but by integrating only over the 
volume between R ns and r = 500 km and adding to this the 
angular momentum that is carried by the gas flux leaving this 
sphere: 

J„s(0 = - ( ° dVpj(r) + A Q JT At' (pjw) | r „J , (2) 

where j is the specific angular momentum and A Q = Anr\. 
This assumes that the asymmetric gas mass outside of r does 
not exert any important torque on the gas mass below r Q . Ac- 
cordingly, we see J ns (f) asymptoting (see Fig.|2]i when accre- 
tion on the NS ends (around 0.6-0.8 s after bounce). At this 
time the asymmetric downdrafts of cool gas filling the vol- 
umes between rising bubbles of high-entropy, neutrino-heated 
matter (Fig. [3] left and right panels) do not reach down be- 
low 500 km, but are replaced by the spherically symmetric 
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Fig. 2. — Time evolution of NS velocity v ns (top-left), angular momentum J ns (bottom-left), latitudinal (top-right) and azimuthal (middle -right) angles of the NS 
kick vector, and 0kkk> respectively, and the angle a s ^ between the NS spin and kick directions (bottom-right). 

neutrino-driven wind around the NS (green region in the right 
panel of Fig. [3). Assuming 7 ns = |J ns | = const after the end of 
our simulations (at t « 1 .4 s p.b.), we obtain a rough estimate 
of the final NS spin period from r sp j n = 27i7 ns /./ ns by consider- 
ing a rigidly rotating, homogeneous sphere of mass M ns , i.e., 
/ ns = |M ns /^ s , with a final radius of R ns = 12 km. 

We find NS kick velocities between 80kms _I (L15-2) and 
405kms~ I (W15-2). In the latter case the acceleration at 1.3 s 
p.b. is still about 150 kms~ 2 (TableQ}. Assuming that half of 
this value applies for another second leads to a final kick of 
nearly 500kms~'. Comparing Figs.[TJand[2]one sees that the 
kick remains very small (< lOkms 1 ) and the kick direction 
fluctuates chaotically in all cases before the explosion takes 
off. Only after the radial expansion of the ejecta sets in and 
the asymmetry pattern of energy and density distribution gets 
frozen in, t he NS acceleration takes place and continues for 1— 
3 seconds (Scheck et al. 2006). Accordingly, the kick finds its 
final direction only shortly afterwards (at 0.4-0.5 s in W15-1 
and W15-2 and 0.6-0.8 s in L15-1 and L15-2). 

Decomposing the radial integral D(9,<p) — J R drp(r) in 

spherical harmonics with the polar axis chosen aligned with 
the kick vector, we find that the I — 1 mode clearly dominates 
after 0.3-0.4 s p.b. in W15-1 and W15-2, the models with the 
largest kick, while this effect is absent in the low-kick models 
L15-1 and L15-2. In Fig. [3] middle panel, the dipolar asym- 
metry of the explosion is clearly visible. The NS recoil (white 
arrow) is opposite to the direction of the fastest shock expan- 
sion (towards the upper left) driven by big high-entropy bub- 
bles. The lower density there is contrasted by much higher 
density and slower gas expansion in the opposite direction. 
Accordingly, the geometrical center of the shock and ejecta 
shell is displaced from the grid center (and NS location) by 
several 1000 km (see middle and right panels of Fig. [3). 

Besides kicks the PNSs attain angular momentum (1- 
2xl0 46 gcm 2 s~ 1 ) without any obvious correlation in magni- 
tude and direction to the recoil. We estimate final NS spin 



periods between 600 ms and 1000 ms (Table [TJ. Evaluating 
D{8, <p) for spherical harmonics with the axis parallel to the 
spin vector, an I — 1, m = 1 mode appears as the dominant 
non-axisymmetric mode in the W15-models after 0.4 s p.b., 
reflecting the mass imbalance on both sides of the spin vector 
as directly visible for W 15-2 in the middle panel of Fig. [3] 

5. PHYSICAL ORIGIN OF THE NS KICKS AND SPINS 

The described hydrodynamic NS acceleration proceeds in 
three stages. (1) Initially convective mass flows and SASI 
sloshing motion of the postshock layer create an anisotropy 
of the mass-energy distribution around the PNS. Convec- 
tive downdrafts, channelling gas accreted through the stalled 
shock into the neutrino-heating region, get deflected to feed 
an asymmetric pattern of high-entropy bubbles. The energy- 
loaded bubbles are created, collapse again, and reappear in 
a quasi-chaotic way to become smaller or larger, absorbing 
less or more neutrino energy. This stochastic bubble forma- 
tion, however, does not cause an appreciable recoil of the NS 
(Fig. |2). (2) When the explosion sets in, the shock and post- 
shock gas begin to expand aspherically, driven by the asym- 
metric inflation of the bubbles. The ejecta gas therefore gains 
radial momentum and its c.o.m. begins to shift away from the 
coordinate origin (Fig. O: The ejecta shell acquires a net lin- 
ear momentum because of different strengths of the explosion 
in different directions. The initial energy and mass asymme- 
try is thus converted to a momentum asymmetry by the con- 
version of thermal to kinetic energy through hydrodynami- 
cal forces. When the expansion timescale becomes shorter 
than the timescale of lateral mixing, the asymmetric ejecta 
structures freeze in. (3) Because of linear momentum con- 
servation, the NS must receive the negative of the total mo- 
mentum of the anisotropically expanding ejecta mass. Hy- 
drodynami c pressure forces al one cannot achieve the NS ac- 
celeration (Sch eck et al.l l2006). As long as accretion down- 
drafts reach the NS, momentum is transferred by asymmet- 
ric gas flows. Stronger accretion on the weaker side of the 
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Fig. 3. — Entropy-isosurfaces of the SN shock and the high-entropy bubbles (left), ray-casting image of the density (middle), and entropy distribution in a 
cross-sectional plane through the center (right) at t = 1.3 s after bounce for model W15-2. The outer boundaries coincide with the shock surface, the viewing 
direction is normal to the plane of the NS kick and spin vectors, which also define the plane for the entropy slice. The kick and spin directions are indicated by 
the white and black arrows, respectively, in the middle figure. The NS location is also marked by a black cross in the right plot and is clearly displaced from 
the geometrical center of the expanding shock. The SN shock has an average radius of 13000 km (a length of 5000 km is given by a yardstick below the middle 
image) but shows a pronounced dipolar deformation, which is clearly visible from the color asymmetry of the postshock gas between the lower right (weaker 
shock with minimum radius of 1 1000 km) and upper left (stronger shock with maximum radius of 15000 km) directions. The middle plot corresponds roughly to 
the projection of the density distribution on the observational plane. Dilute bubble regions are light-colored in white and yellow, while dense clumps appear more 
intense in reddish and bluish hues. The blue circle around the NS represents the dense inner region of the spherically symmetric neutrino-driven wind. This wind 
is visible in green in the right image and is bounded by the aspherical wind termination shock. The wind is shocked to higher entropies on the left side, where it 
passes the termination shock at larger radii because of the faster expansion of the preceding SN ejecta. 

blast and more mass loss in the neutrino-driven wind on the 
other side cause a recoil opposite to the main explosion di- 
rection. However, the largest (time-integrated) kick contribu- 
tion, which continues even after accretion has ended and the 
wind has become spherical (Fig. |2j, results fr om the gravi- 
tation al pull of the anisotropic shell of ejecta (IScheck et alJ 
2006). A hemispheric asymmetry of the mass distribution of 
Am = +10" 3 M o in a shell expanding away from the NS from 
a radius r; = 100 km with v s = 1000 km s _1 can drag it to a 
velocity of v ns « 2GAm/(r;v s ) a 2700 km s -1 . Ejecta asym- 
metries can thus effectively mediate a long-lasting pull on the 
NS. 

We do not expect any significant contribution to the NS 
kick by anisotropic neutrino emission, because most of the 
neutri no energy is radiated from the spherical neutrino spheric 
layer (Sch eck et al.l|2006l) . Moreover, including the NS mo- 
tion is unlikely to change our results, at least statistically for a 
larger sample of models. This was concluded by lScheck et al.1 
(2006) from tests in which the gas surrounding the NS in 
the grid center was allowed to move with (— v ra ) by apply- 
ing a Galilei transformation. The displacement of the NS 
is small compared to the asymmetry scale of the ejecta that 
cause a grav i tation al acceleration of the NS over seconds. 
IScheck et aD ((2006) also observed little influence of the NS 
motion on the neutrino emission asymmetry and its a ssoci- 
ated (v ery small) NS recoil, in contrast to findings by iFrverl 
(12001 . 

Our 3D simulations were followed much longer b eyond 
the on set o f the explosion than pr evious 3D studies by IFrverl 
(2004) and iFryer & Youngl (120071) . This is essential because 
the hydrodynamical NS ki cks are no impulsive p re-explosion 
effect (as, e.g., assumed by Janka & Miiller 1994) but a long- 
time post-explosion phenomenon. The global asymmetry of 
the mass-energy distribution in the postshock gas created be- 
fore the explosion mainly by SASI activity (as analysed in 



detail by IScheck et"ail (J2008) for the considered SN core 
conditions) is converted to a momentum asymmetry of the 
ejecta and the NS only gradually after the initiation of the 
blast. Since the acceleration continues for seconds, kicks of 
>500kms _1 are possible although an extreme dipolar shock 
deformation and persiste nt one-sided accretion as seen in 
some 2D models of lScheck et al.l (|2006) are absent in the pre- 
sented 3D results. 

The NS spin grows considerably faster than the kick but 
saturates also earlier (Fig. |2j. Angular momentum is trans- 
ferred to the NS by off-center impacts of the anisotropic 
accre tion downdrafts dBurrows et al.lfT995t ISpruit & Phinneyl 
1998). When the accretion ends, the exterior gas separates 
from the compact remnant. Different from mediating an ac- 
celeration by their gravitational drag, the anisotropic ejecta 
then cannot exert any torque on a spherical, point-like com- 
pact remnant. Therefore the NS spins remain fairly slow, be- 
cause in spite of an t — 1 , m = 1 density asymmetry we do not 
find the strong, ord ered rotational gas flows near the accret- 
ing PNS as seen by Blondin & Mezzacappa (2007). Possibly 
the investigated SN core conditions are not favorable for the 
growth of rotating spiral modes (see, e.g., Fernandez 2010) or 
their growth timescale is longer than the explosion timescales 
in our models. 

6. CONCLUSIONS 

Our 3D results support th e NS recoil scenario propo sed on 
the basis of 2D models by IScheck et all d2004l 2006). The 
still small set of 3D simulations for two 15M stars yields 
kicks between ~50km s _1 and ~500km s _1 in the range of the 
measured space velocities of young pulsars. Since even for 
our most extreme model the explosion asymmetry of the SN 
is still fairly small, NS acceleration well beyond 500 km s _1 
seems possible. But since the kick mechanism is a quasi- 
chaotic and stochastic phenomenon, rare cases with more than 
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1000 km s as in 2D (lSchecketal.1 12006) will require many 
more model runs, varying the initial random seed perturba- 
tions as well as the stellar progenitor and explosion energy. 

This will also be necessary for statistically significant con- 
clusions on the relative orientations of NS spins and kicks. 
Considering nonrotating stars we find spin periods of -500- 
1000 ms, but do not see reasons for a spin-kick alignment. 
For a closer assessment of this question, however, also the 
investigation of progenitors with rotation is necessary. Core 
rotation due to the angular momentum of the progenitor (or 
established by a first strong thrust) might impose a prede- 
fined direction and lead to rotational averaging of the fol- 
lowing NS acceleration so that a preferentially aligned spin- 
kick distribution emerges instead of random distributions of 
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away from the NS motion. 
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